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Abstract 



■ The Dickman function F{a) gives the asymptotic probabihty that a large integer N has 

, no prime divisor exceeding N"'. It is given by a finite sum of generahzed polylogarithms 
defined by the exquisite recursion Lk{a) = — jj dxLk_i{x/{l — x))/x with Lo(a) = 1. 

(— I ' The behaviour of these Dickman polylogarithms as a — > defines an intriguing series of 

Ci-i' constants, C^. I conjecture that exp(7z) /r(l — z) is the generating function for Yl,k>Q ^k^^ ■ 

r~| ! I obtain high-precision evaluations of F{l/k), for integers k < 11, and compare the 

' Dickman problem with problems in condensed matter physics and quantum field theory. 

1 Introduction 

> 

The fraction of positive integers less than with no prime divisor greater than N°' tends 
"/^ ■ to a finite limit F{a) as — > oo. Clearly F{a) = 1 for a > 1. For < a < 1, the 

Dickman [21] function F{a) satisfies the remarkable differential equation [18\ [35l [TBI [26] 

o: ^^"^"a^vr-a. 



F'{a) = If (-^ ] (1) 



and may be computed as a finite sum of terms 



K{a) 

F{a) = Lu{a) (2) 



; fc=o 

where Lq^o) = 1, K{a) is the largest integer such that aK{a) < 1, and the recursion 



Lk{t) = -r Lk-i { j^) ^ (3) 



defines the Dickman polylogarithm of weight A; > as an iterated integral. I shall prove 
two theorems that allow me to obtain at least 100 digits of Lk{t) for all weights A; < 10 
and hence shall conjecture the generating function for the Dickman constants 

C, = ynni^LM-i^i^) (4) 

where Cq = 1 and is determined by the boundary condition ^fc(T) = for A; > 0. 
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2 Dickman trilogs from standard trilogs 



Since Li{t) = log(t), we have Ci = and F(l/2) = 1 — log(2). Moreover, we easily 
obtain [17] 

_ze! 

12 



where Lifc(t) = J2n>o^"^ /''^ t'^^ standard polylogarithm of weight k and C2 = 
ensures that L2{^) = 0. Thus 

°° 1 1 7r2 

F(l/3) = l-log(3) + E;^ + 2l°g'(3)-Y^ (6) 

n=l 

pa 0.048608388291131566907183039343407421354329580478141 (7) 
may be computed with great facihty. 

After considerable effort, I shall prove in Section 5 that 

L3(t) = Li3(2t-l)-Li3(l-2t)-Li3(t)-Li3(2-l/0 

+ (1^) ^^^^^ + T " I + i 

and hence obtain the Dickman constant 

C3 = lim (^Lsit) + ^ log(t) - ^ log3(t)^ = -i C(3). (9) 

It follows from (j2l5|8p that F{a) reduces to standard polylogs for a > ^. In particular, 
the asymptotic probability that an integer has no prime divisor greater than A 4 is 

F(l/4) = l-21og(2)+Li2Q)+21og2(2)-^ 

- Li3 (^) - Li2 log(2) - ^ log3(2) + g C(3) (10) 
^ 0.0049109256477608323527391509236151860324842974176929 (11) 
with 800 good digits obtainable in less than 10 milliseconds. 



3 Conjecture for the Dickman constants 

I shall show, in later sections, how to evaluate F{a), for j > a > jq, and the intriguing 
constants C^, for 3 < k < 10, by one-dimensional numerical quadrature, at precisions well 
in excess of 100 decimal digits. This is amply sufficient to obtain very reliable conjectures 
for the analytic form of C^. The complexity of these calculations is in marked contrast to 
the final simplicity of the conjectured results. For example, 

_ 1 u. a) - ^ Li3 f i) + Lb a) + ?1M5M 

8 V4/ 4 V4/ 12 V4/ 8 

7r2log^(2) ^2i^g(2)log(3) log=^(2)log(3) 51og^(2) 

24 6 2 8 ^ ' 
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has a very simple conjectural evaluation as j^Cf, which holds at 800-digit precision. The 
complexity of evaluation (jl2p results from delicate integration by parts, so as to remove 
divergent terms at the lower limit of integration as t ^ in ([3]) and leave a finite integral 
whose integrand contains no trilogs. 

The constant Cg is of particular interest, since that is the first case where a multiple 
zeta value (MZV) |40 t [9l[T0 t 11 ^ 16] might have occurred. However, after far more intricate 
numerical quadratures, the LLL [29j algorithm implemented by the Pari-GP [19j procedure 
lindep gave the conjectured evaluation 

^ C(5)C(3) _ [C(3)vr]^ _ 67^8 
^ 15 216 29030400 ^ ^ 

which is free of the irreducible MZV of lowest weight, namely C(5, 3) = J2m>n>o ^/im^n^)- 
After noticing that —67 is the first non-unit integer in the sequence 1, —1, 1, —1, —67, — 1, 
given by Neil Sloane p7j in OEIS entry A008991, for numerators of coefficients in the 
expansion of 



/ sin(x) _^ 67x^ ^L_^o(x^A 

V X 12 1440 24192 29030400 5677056 V / ' ^ ' 

I was led to conjecture the wonderfully compact generating function 

G{z) ^ ^^pM 1 V Cz'^ (15) 

r(i-z) ' 

where 7 is Euler's constant. The equivalent, yet more illuminating, formula 



neatly accounts for Sloane's tell-tale integer 67 in ()13p and leads to the conjecture 

^ _C(9) vr^ _ vr^ _ mf .^7. 

^ 9 84 7200 72576 162 ^ ' 

which was tested by numerical quadrature of a product of Dickman tetralogarithms, 
thanks to the following theorem. 



4 Integration by parts 

I define an auxiliary family of functions by the recursion 

, N fn(x) dx , , 

/n+1 t = - / ^ , 18 

Jt x\ \ — nx) 

with /o(t) = 1. Then f\{t) = log(t) = Li{t) and the dilogarithm 
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is easily related to L2{t), using integration by parts. 
Theorem 1: Let 



. ^ / N /" ^ / X \ fn(x) dx , , 

M.,,(,).-/ W.( i_(„_,i)J ^I^ (20) 
for any pair of integers k > n > 0. Then this integral evaluates to 



m=0 



Proof: At n = 0, definition (|2Up gives Mi.fi{t) = Lk{t)^ by virtue of recursion Q for 
Lk{t). Hence (pT]) holds at n = 0. For A; > n + 1 > 0, recursion ([TS]) for fn+i{t) allows an 
integration by parts in definition (|20p . to obtain 

Mk,n{t) = Lk-n-l (- T^—TT:) fn+l{t) " Mfc,„+i(t) (22) 

VI - (n + l)t/ 

with a vanishing constant term, since the Lfc_„_i term vanishes at t = 1/k, where its 
argument t/(l — (n + l)t) evaluates to l/(/c — n — 1). Hence I prove (12ip by induction, for 
all /c > n > 0. 

Comment: This is a very powerful result, peculiar to the Dickman problem. For 
example, it allows us to compute the Dickman heptalogarithm Li{t) very accurately, as 
a single integral of a product of trilogarithms, instead of having to evaluate a four-fold 
iterated integral of trilogs. 

The key to the method is to observe that fn+i{t) = M„+i^„(t). To prove this, I set 
/c = n + 1 in the definition ([20]) of Mk^n{t) and then use recursion ([TH]) for fn+i{t). I thus 
prove the claimed result ()19p for /2(t) by setting k = 2 and n = 1 in the evaluation ()21l) 
of Theorem 1. Similarly, yet much more importantly, I obtain 

hit) = log (y^) ^2^*) " ^2 (r^) + ^'^^^ ^^^^ 

by setting k = 3 and n = 2, thereby avoiding duplication of the considerable effort 
expended in obtaining the trilogarthmic result ([8]) for L^it). Then, by setting k = 7 and 
n = 3 in the theorem, I obtain the Dickman heptalogarithm 

LAt) = i«(r^)M<)-i5(^)AW + L.(^)/3W 



+ /'i3(:r^|#i^ (24) 



t \1 — Ax J x{l — 3x) 

with a one-dimensional quadrature of products of known trilogs, in the final term. More- 
over, there are two methods of evaluating LQ{t) as an integral of products of known dilogs 
and trilogs, using the pair {L2, fa) or the pair (^3,72) in the final integrand. For L5(t) 
there are three methods, using (L3,/i), {L2, f2) or {Li,/^). Of these, the (-^25/2) pair 
is the most efficient. For Li{t) there are four methods, using (-L3,/o), {L2,fi), (Li,/2) 
or (Lo,/3), where Lo(*) = /o(*) = 1- Efficiency dictates that one should use either the 
second or third pair; caution suggests that one should use all four methods to check the 
accuracy of numerical quadrature. 
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5 Dickman heptalogs from standard trilogs 



To justify this methodology, I must first prove the claimed result ([8]), which reduces the 
Dickman trilog L3 to standard trilogs and establishes my claim that C3 = — |C(3). I begin 
by proving the rather simple identity 

L3{t) = C3 + C2 log(t) + I log3(t) - Li3(t) + Li2(t) log(t) 

2 

where C2 = —■j2 determined by ([5]) and C3 is an integration constant, to be deter- 
mined later by the requirement that -^^3(^) = 0. The proof of (|25p is symptomatic: we 
simply differentiate with respect to t and check that L'^{t) = L2{t/{\ — t))/t, as required 
by ([3]). Here, as ever, I use the relation Li'^(t) = Lifc_i(t)/t, with Lii(t) = — log(l — t). 

Next comes a harder part, namely to perform the integration in (j25p . Using integra- 
tion by parts, I was able to reduce the problem to an instance of equation 8.111 in the 
fascinating and enormously informative book by Leonard Lewin [30j . long since sadly out 
of print. Here, I am content to state and then to prove that 

+ Li3(l - 2t) - Li3(2t - 1) + Li3 (3-3^) 

+ Li2(t) + ^ Li2(2t) - ^ Lii(2t) Lii(l -t) + ^ Li2(l)) Lii(2t) . (26) 

This is clearly true at t = 0, since Li3(l) = C(3) and Li3(— 1) = — 1(^(3). Thus it is sufficient 
to show that the right hand side of (j26p has a vanishing derivative. This derivative is of 
the form Di{t)/t + D2{t)/{1 — 2t), where Di and D2 are rather complicated combinations 
of dilogs and products of logs. However, it is easy to show that -Di(O) = D2{0) = 0. 
Hence, to prove (j26|) . it is sufficient to show that D[{t) = D'2{t) = 0, which may be done 
by elementary manipulation of logs. 

Inverting the argument of Li3(— t/(l — 2t)), in ()26p . and imposing the boundary con- 
dition L3(i) = 0, in ([25]), I arrive at the claimed result ([8]) for L3{t) and the claimed 
evaluation 3C3 = — C(3), provided that 

3 (2Li3 Q) - Li3(-3)) - log=^(3) = y C(3). (27) 

The final hurdle of proving (j27p was the most challenging. Spencer Bloch and Herbert 
Gangl told me that they expected the combination 2Li3(|) — Li3(— 3) to evaluate to some 
rational multiple of C(3), "modulo logarithms". Yet it appears that such rational numbers 
are as distressingly hard to derive from first principles as they are disturbingly easy to 
guess from low precision numerical computation. My claim that C^, = — ^C(3) requires 
this rational multiple to be with a denominator divisible by 3. It was rather hard to 
see how a simple functional equation for trilogs of a single variable might produce such 
a denominator. Accordingly, I resorted to the Spence-Kummer functional relation for 9 
trilogs of a pair of variables, given in equation 6.107 of Lewin's book [30]. Setting x = — 1 
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and y = I in that ornate identity and inverting the argument of Li3(— |), I proved (p7|) . 
obtaining the combination 2Li3(l) — 6Li3(— 1) on the right hand side. Thus the "morahy 
rational" coefficient of C(3) for the combination 2Li3(^) — Li3(— 3) is now proven to be 
|(2 — 6(— |)) = thereby rescuing this particular problem from what Sasha Beilinson [3] 
memorably referred to as "the burdock thicket of generalities" . 

Having thus proven the trilogarthmic input for the method of Theorem 1, I am able 
to compute F{a) for a > |, with great ease, and have provided 800 good digits for 

F(l/5) ^ 3.5472470045603972983389451077062356095164361057262 x lO""^ (28) 

F(l/6) ^ 1.9649696353955289651754986129204522894596719809623 x 10"^ (29) 

F{l/7) f« 8.7456699532939166955802835727699721733804719764580 x 10"^ (30) 

F(l/8) « 3.2320693042261037725997853617282161576194751628024 x 10~® (31) 



at the web page http : //physics . open. ac .uk/f "dbroadhu/cert/smoctic . txt 



A significant amount of LLL analysis produced no simple integer relation of F{l/5) 
to values of standards polylogs and their products. It might be interesting to investigate 
this issue more intensively, using David Bailey's parallelization [3] of Helaman Fergu- 
son's PSLQ [22] algorithm, since more than 800 digits of -L4(l/5) may now be computed 
from (l3|), using the explicit trilogarithms in ([8]). My own opinion, however, is that the 
Dickman probability F(l/5) does not reduce to standard polylogs. 



6 Dickman octalogs from Dickman tetralogs 

It was rather frustrating that Theorem 1 did not provide a method for computing Dickman 
octalogs as single integrals of products of standard polylogs, since it was precisely at weight 
8 that I wished to have high-quality numerical data with which to determine whether the 
first MZV, namely C(5, 3), might show up in Cg. It was moral support from Mike Oakes |34j 
that determined me to push the investigation above weight 7. 

The barrier that must now be surmounted stems from the fact that 

^^^f"' ^ i" ('°^ (its) ■''^(^■' + 5 '°^'<-' "'(-^-)) WTT) 

appears in the Dickman tetralog, since L^{t) — C4 + M^{t/{1 — 2t)) may be reduced to 
standard polylogs and their products. It is an easy matter to evaluate a single instance 
of ()32p to high precision. For example, MA(h) in (pi) was computed to 800 digits, in order 
to determine the value of C4 that results from the boundary condition L4{j) = 0. However, 
it is quite another matter to obtain good results for Cg and Cg, which are integrals with 
M4 in their integrands. When t is significantly less than |, one may efficiently evaluate 
M/j^lt/ (1 — 2t)) by using Taylor series for log(l + 2x) and 1/(1 + x) under the integral sign 
in ([52]) . Yet to compute Cg and Cg we need Li{t) for all arguments < t < j, inside 
the integrals for Lg and Lg. So the remaining problem involves an investigation of the 
behaviour of L4 (t) in the neighbourhood of t = ^. 

Theorem 2: Let gniz) = Ln{^ — z) for n > 0. Then gn{z) has a Taylor series with 
rational coefficients, beginning with gn{z) = GnZ"^ + 0{z^'^^), where Gn = (— l)"'(n"'/n!)^. 
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Proof: We know that gi{z) = log(l — z) = — z + 0(z^) has a Taylor series with rational 
coefficients. Setting k = n + 1 in recursion ^ and transforming variables, I obtain 



. ^ r {n + l)h \ jn + l)dy 

gn+i[z) = - 9n\^-^ — 7 — — p-h^ 7 — r-p-- (33 

Jo + n{n + l)y J 1 - (n + l)y 

Now suppose that the claim is true for gn- Then, by binomial expansion under the 
integral sign in (p3]) . it is also true for gn+i, since Gn+i/Gn = —{{n + l)^/n^)". Hence, 
by induction, the claim is true for all n > 0. 

Comment: Starting with gi{z) = log(l — z), I used three iterations of recursion (I33p to 
compute sufficient terms in the rational Taylor series for 54(2) = ^z^ + 0{z^) to achieve 
at least 240 digits of precision for L^it) = g^^j — t) in the region ^ > t > = 0.2358. 
Then, for t < ^4, I proceeded as follows. 



Let a„ = {-bn + (-l)"c„)/n, where 



k=l 



and Cn is the coefficient of z" in the Taylor series for log(l + 2z)Li2(2;)/(l + z). Then, 
with x = t/{l — t) and y = t/[l — 2t), the summation in 

00 

L4(t)-C4 = ^{a„ + 5„log(y)}(-yr + L3(x)log(t)+Li2(-?/)log2(y) 



n=2 

Li2(-x) + \ log2(x)) (^^ - i log\y)^ + ^ log^y) (35) 

provides a viable method for computing L4^{t) for t <ti = 0.2358. 

A sanity check was provided by comparing the two strategies at t = ^4 and verifying 
that 240 good digits of the conjectured value were obtained for C4. To achieve this, 
I expanded 54(2) up to 0{z^^^) and took 4800 terms in the summation over powers of 
—y = —t/{l — 2t) in (j35p . The relatively small number of terms used for g^ resulted from 
the fact that it is rather laborious to perform three rational binomial iterations of recur- 
sion (j33|) . By storing coefficients of the expansions, I am able to evaluate 240 good digits 
of the Dickman tetralog L4 in less than 50 milliseconds and hence can efficiently com- 
pute the Dickman octalog, Lg, and the nonalog, Lg, as one-dimensional quadratures with 
integrands that call the procedure for L4. By these means, I arrived at conjecture (jl3p 
for Cs and then, thanks to OEIS sequence A008991, inferred the wonderfully compact 
conjectured generating function exp(72;)/r(l — z) for J2k>o^kZ^ 1 which then gave con- 
jecture (fT7|) for Cg, also now verified at high precision. Finally, I record 50 good digits of 
the Dickman probabilities 

F(l/9) ^ 1.0162482827378365465348539356956957838244399586581 x 10"^ (36) 
F(l/10) ^ 2.7701718377259589887581212006343423263430066501156 x 10^1^(37) 

noting that 6 good digits were given in [32], which corrected serious errors in [5]. 
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7 Context and conclusion 



Thus far, I have spared the gentle reader exphcit reference to quantum field theory (QFT) 
or to condensed matter physics. I trust that s/he will now permit me to reveal the physics 
context for investigating Dickman's mathematical |2H [THl [35} [161 1271 l26l [H [12] problem. 

The collision of a pair of protons in the large hadron collider (LHC) [31], at energies 
never before achieved in a particle accelerator, is described by QFT in terms of the collision 
of quarks, or gluons, carrying a fraction x of the momentum of a proton. The key thing 
that we need to know, to make sense of the possible outcomes of these LHC collisions, 
is the probability distribution function (PDF) in x. This is not known a priori; rather 
it is inferred from electron-proton collisions at lower energies, where half of the problem, 
namely the electron, is already well understood. Yet the input PDF, from electron- 
proton collisions, is not immediately usable at the LHC. Rather, it must be "evolved", 
up to the higher LHC energy. The procedure for doing this is completely understood, 
in principle, yet involves formidable mathematical challenges, in practice [21 [39]. It is 
performed as a perturbation expansion in the coupling constant of the strong interaction 
between quarks and gluons. At successive orders in this compelling expansion, generalized 
polylogarithms [361 [24] in x make their appearance. 

From the perspective of QFT, the Dickman problem might seem to be rather routine: 
as soon as the polylogarithmic structure of the recursion ([3]) for is exposed, it is clear 
that the formidable technical machinery of QFT ^[24l[25l[33l[6l[Ml[20l[28jhas much to 
offer to the computation of this problem in number theory. The fascinating circumstance 
is that each iteration ([3|) of the Dickman polylog brings its own distinctive constant 
Ck in ©• This also mimics QFT, where limiting values of generalized polylogarithms are 
likewise zeta values, until one reaches weight 8, when a multiple zeta value appears [141 
[T3l [151 [6] . The Dickman problem would have shared even more features with QFT had Cg 
contained this first irreducible MZV, namely C(5, 3) = J2m>n>o ^/(m^n^)- In this paper I 
have shown, with overwhelming probability, that ^(5, 3) does not appear in Cg. Moreover, 
I have conjectured that is, most wonderfully, the coefficient of in the expansion of 
exp(7z)/r(l — z) and hence zeta- valued for all k. 

This does not, however, mean that the Dickman problem now lacks interest for physi- 
cists. On the contrary, it resembles fascinating problems in condensed matter physics. 
In the study of quantum spin chains ff] , hugely demanding calculations gave zeta- valued 
results up to weight 7. When Dirk Kreimer and I asked Valdimir Korepin to climb the 
next cliff, up to weight 8, we genuinely did not know whether an MZV would appear. In 
fact it did not [8j. Thus the Dickman problem sits very well with the study of spin chains. 

In conclusion, I suggest that the current investigation has added to our understanding 
of the Dickman function and has also confirmed how much more demanding is the poly- 
logarithmic structure of QFT, which appears to surpass the Dickman function of number 
theory, and the spin chains of condensed matter physics, both in its grand challenge and, 
I believe, in its eventually to be comprehended great beauty [23] . 
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